Stable nontrivial topology in ultrathin Bi (111) films: a first-principles study 
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Recently, there have been intense efforts in searching for new topological insulator (TI) materials. 
Based on first-principles calculations, we find that all the ultrathin Bi (111) films are characterized 
by a nontrivial Z2 number independent of the film thickness, without the odd-even oscillation of 
topological triviality as commonly perceived. The stable nontrivial Z2 topology is retained by the 
concurrent band gap inversions at multiple time-reversal-invariant k-points and associated with the 
intermediate inter-bilayer coupling of the multi-bilayer Bi film. Our calculations further indicate 
that the presence of metallic surface states in thick Bi(lll) films can be effectively removed by 
surface adsorption. 

PACS numbers: 73.21.Ac, 71.70.-d, 73.20.At 
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As a new insulating phase in condensed matter, the 
topological insulator (TI) has recently attracted a great 
deal of attention [l|. The TI is distinguished from the 
conventional insulator by its unique gapless surface states 
residing in the middle of band gap as a consequence of 
the so-called Z2 topology encoded in the wavefunctions. 
There has been an intensive search for TI materials. Al- 
though quite a few compounds have been found to be 
3D TIs[2-l6|, up to now only the HgTe quantum well is 
verified to be a 2D TI (or a quantum spin Hall insulator) 
experimentally Recently, Murakami predicted a sin- 
gle bilayer (BL) Bi (111) film to be an elemental 2D TI 
[1, and further speculated the multi-BL Bi (111) film 
to exhibit an odd/even oscillation of topological triviality 
with film thickness 3 by considering the multi-BL film as 
a stack of BLs with no or very weak inter-BL coupling 
(loj . Generally speaking, the 2D TI phase in thin films is 
commonly perceived to depend sensitively on film thick- 
ness, as shown for the Bi2Se3 and Bi2Te^ ultrathin films 
lll | . The requirement of fine tuning the thickness makes 
the experimental fabrication of 2D TIs rather difficult. 
Therefore, it is desirable to search for new materials or 
new schemes to obtain a 2D TI. 

In this Letter, based on the first-principles calculations 
of band structures and wavefunction parities, we find that 
in fact all the Bi (111) ultrathin films are characterized 
by a nontrivial Z2 number independent of the film thick- 
ness. The films with 1- to 4-BL are intrinsic 2D TIs, 
while those with 5- to 8-BL are 2D TIs sandwiched with 
trivial metallic surfaces that can be extrinsically removed 
by surface adsorption. This finding is in direct contrast 
to the odd/even oscillation of topological triviality spec- 
ulated for the Bi films Q, as well as to the thickness- 
dependent topology shown for other 2D films 11|. The 



surprisingly stable 2D TI phase in the Bi(lll) films are 



found to be retained by a concurrent band gap inver- 
sion at multiple time-reversal-invariant k-points (TRIKs) 
when the film thickness is increased. A detailed analysis 
of the 2- and 3-BL Bi films indicates that the intermedi- 
ate inter-BL coupling plays an important role in defining 
their unique topological property. 

Bismuth is one of the main group elements that has 
the strongest spin-orbit coupling (SOC), a fundamental 
mechanism to induce the Z2 topology. For this reason, 
many 3D TIs make use of Bi, even though 3D bulk Bi 
itself is topologically trivial [2| . The electronic properties 
of Bi, such as the bulk band structure 12|, the surface 
state s I13II and the semimetal-to-semiconductor transi- 
tion [IJ, |l5| , have been well established in the literature. 
The outermost shell of Bi has the electron configuration 
6s^6p^. It tends to form three bonds to close the shell. 
The single-crystal Bi has a bi-layered structure, with an 
ABC stacking sequence along the (111) direction (Fig. 
la). Within each BL, every Bi atom forms three a bonds 
with its nearest neighbors in a trigonal pyramidal ge- 
ometry. Projecting onto the (111) plane, the BL forms 
a hexagonal lattice with two atoms per unit cell (Fig. 
lb). There are three key structural parameters to define 
the lattice: the in-plane lattice constant a, the intra- 
BL bond angle a (or the intra-BL height difference di), 
and the inter-BL spacing ^2. Our calculated Bi crystal 
structural parameters as shown in Fig. la and b agree 
well with previous calculations [l3|. The ultra-thin Bi 
(111) thin films, consisting of a few number of stacked 
BLs, have slightly relaxed structural parameters relative 
to the bulk values. Specifically, a approaches 90° and 6,2 
increases by about 6%. This implies a slightly enhanced 
p-orbital feature of the intra-BL bonds and a weakening 
of the inter-BL coupling. 

To identify the 2D TI phase, a single Z2 topological 
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FIG. 1: (a) The hexagonal unit cell of single-crystal bismuth, 
(b) The top view of the Bi lattice, (c) The first Brillouin zone 
of the hexagonal lattice. 

number (v) is used as the "order parmeter" v = 1 

indicates a topologically nontrivial phase; v — Q indicates 
a trivial phase. The calculation of the Z2 number can be 
dramatically simplified by the so-called "parity method" 
0, if the system is space inversion invariant, as the case 
of the Bi (111) film. Accordingly, the Z2 number of Bi 
films can be obtained from the wavefunction parities at 
four TRIKs {Ki), one F and three Ms, as 



TABLE I: The total parity at the F and M points and the Z2 
number of Bi (111) films with different thickness 
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where ^ = ± is the parity eigenvalues and iV is the num- 
ber of the occupied bands. 

Single-electron wavefunction parities are calculated 
within the density functional theory using the plane wave 
basis, as implemented in the ABINIT package II6II . We 
employ the local density approximation (LDA) [l7[ and 
the Hartwigsen-Goedecker-Hutter pseudopotential [l8| . 
which is generated on the basis of a fully relativistic all- 
electron calculation and tested to be accurate for heavy 
elements like Bi. The spin-orbit coupling is included in 
the self-consistent calculations as described in [l^. To 
model the thin film, a supercell of slab is used with peri- 
odic boundary conditions in all three dimensions with a 
10 A thick vacuum layer in the (111) direction to elim- 
inate the inter-slab interaction. A plane wave cutoff of 
24 Ry and a F-centered k-point mesh of 10 x 10 x 1 are 
used. All the atomic positions are fully relaxed for each 
film. 

Table 1 shows the calculated numbers for 1- to 8- 
BL films. Surprisingly, all the films we calculated are 
characterized by the nontrivial Zi number {y = 1), in 
direct contrast to the oscillation of topological triviality 
as commonly perceived. We notice that the total par- 
ity (5(r) and (5(M) simultaneously change their signs for 
every 2 BLs. This kind of parity oscillation, which orig- 



FIG. 2: (a)The direct (F-point) and indirect band gap as a 
function of the film thickness, (b) The band structure and 
the parity information of 5-BL. (c) A schematic depiction of 
the electronic property of an Bi (111) film before and after 
the H-termination. (d) The band structure and the parity 
information of H-terminated 5-BL. 

inatcs from the inverted band gap of a 3D TI, has been 
reported in Bi2Se3 and Bi2Te3 ultrathin films, but only 
at the F point jT2|. The uniqueness of Bi film is the in- 
verted band gaps both at the F and M points due to the 
strong SOC of Bi Q. As a consequence, the parity oscilla- 
tion under the quantum confinement occurs concurrently 
at all TRIKs (one F and three Ms) as the film thickness 
increases. Being the product of S{T) and 5{M), the Z2 
number shows up as the "beat" of two oscillations, which 
makes the 2D TI phase in Bi ultrathin films much more 
stable than in the Bi2Sez and Bi2Tez films. We expect 
that the Z2 number of Bi films will eventually change at 
some point when the phase difference between the oscil- 
lations at F and M accumulates to tt, but this requires 
calculations of much thicker Bi films possibly beyond the 
current computational capability. 

In Fig. 2a, we plot both the direct and indirect band 
gaps as a function of the film thickness. Within the calcu- 
lated film thickness range, the direct gap always remains 
open, which is essential for a well-defined Z2 number. 
Below 4 BLs, the film is a semiconductor having a non- 
trivial Z2 number, and hence representing an intrinsic 
2D TI. However, the indirect band gap becomes nega- 
tive above 4 BLs, leading to a semiconductor-semimetal 
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transition[15j. The semi-metallization arises from two 
overlapping bands around the Fermi level, as shown in 
Fig. 2b for the 5-BL film as an example. It has been 
pointed out 13|, |l5| that these two bands have evident 
surface band features and tend to become gapless at the F 
and M points in the limit of a semi-infinite system. From 
this view, the films from 5- to 8-BL can be regarded as a 
2D TI sandwiched between two trivial metallic surfaces 
(top panel, Fig. 2c). The meaning of "trivial" here is 
two-fold. On the one hand, if we consider the surface as 
an individual 2D system, its number is (trivial), as 
obtained from the surface band parities at the F and M 
points (see Fig. 2b). Therefore, the surface bands have 
no contribution to the nontrivial number of the film. 
On the other hand, in the limit of semi-infinite system, 
because bulk Bi is a 3D Zi topologically trivial insula- 
tor, these trivial metallic surface bands are not robust 
as those of a 3D strong TI and hence, in principle, can 
be easily removed by surface defects or impurities, e.g. 
via surface adsorption. To test this idea, we have ter- 
minated the two surfaces with H atoms as schematically 
depicted in the bottom panel of Fig. 2c and repeated 
the calculation. We find that upon surface adsorption, 
the two surface bands are separated apart opening a gap 
around the Fermi level, as shown in Fig. 2d. Meanwhile, 
the H atoms introduce additional occupied bulk bands, 
which are found to be topologically trivial, so that the 
Z2 number remains nontrivial (see Fig. 2d). Thus, by 
the extrinsic effects of H surface adsorption, the thicker 
films above 5 BLs are effectively converted into 2D TIs 
similar to those ultrathin films below 4 BLs. We note 
that because LDA is known to underestimate the band 
gap, the actual transition thickness is likely to be thicker 
than the 4-to-5 BLs, but the overall trend of behavior we 
show here should remain valid. 

It has been predicted that the 1-BL Bi (111) film is 
a 2D TI [i, i]. If we could regard the n-BL film as a 
stack of these nontrivial 1-BL films and without inter- 
BL coupling, all the bands would have n-fold degeneracy 
and there would be naturally an odd-even oscillation of 
Z2 topology: V = \ loY the odd-BL stacks and v — ^ 
for the even-BL stacks [l3| as every additional BL fiips 
the Zi number. Such odd-even oscillation of Z2 topology 
from the zero inter-BL coupling limit can be extended to 
the weak inter-BL coupling regime under the adiabatic 
approximation, which was speculated to be applicable to 
the Bi(lll) films fH. However, our first-principles results 
show that the Bi (111) film may not be adiabatically 
connected to the zero inter-BL coupling limit. 

The Bi (111) film actually represents a special inter- 
esting class of films having an intermediate inter-BL cou- 
pling strength. The inter-BL bond energy is calculated 
to be 0.3-0.5 eV per bond, which is noticeably larger 
than the typical values of weak interfacial bonds, such 
as Van der Waals bond, but smaller than the values of 
typical chemical bonds. This intermediate inter-BL cou- 
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FIG. 3: The band structure of (a) 2-BL and (b) 3-BL under 
different inter-BL spacing, (c) The energy level at the F point 
of (left) 2-BL and (right) 3-BL as a function of the inter-BL 
spacing. The red (circle) lines indicate even-parity levels. The 
blue (square) lines indicate odd-parity level. 



pling may have a significant influence on the topological 
property. To reveal the influence, we have calculated 
a set of "model" Bi films with their inter-BL coupling 
tuned gradually from the real intermediate regime to the 
hypothetical weak coupling regime that could be adia- 
batically connected to the zero coupling limit. This is 
done by artificially increasing the inter-BL distance from 
the equilibrium value c?2 with an increment of Ad2 7 using 
the supercell technique. Here we take 2-BL (3-BL) film 
as an example of even (odd) number of BL films. Fig- 
ure 2a and 2b show the evolution of band structures as a 
function of Arf2 for the 2-BL and 3-BL films, respectively. 
For A(i2 — 0, the figures show the realistic band struc- 
tures of the 2-BL and 3-BL films, while for A(i2 = 4A, 
the corresponding band dispersion of hypothetical films 
shows a weak inter-BL coupling case that is adiabatically 
connected to the zero coupling limit, based on topologi- 
cal analysis. We note that as the A(i2 is tuned from 4A 
to OA, the direct gap at F point undergoes a closing and 
re-opening process, indicating that the realistic band dis- 
persion of Bi films may not be adiabatically connected to 
the zero coupling limit. 

To trace the change of topology with the inter-BL cou- 
pling, we analyze the parity property of the energy levels 
at F point as a function A(i2 in Fig. 2c for the 2-BL 
(left) and 3-BL (right) Bi films, respectively. When Ac?2 
is reduced, the n-fold degeneracy for the zero coupling 
limit as mentioned above is lifted by the inter-BL cou- 
pling. Although the variation of band gaps is similar 
for the 2-BL and 3-BL films, the change of total parity 
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of all the occupied bands shows quite different behav- 
iors. For the 2-BL film, each doubly degenerate level 
is split into one odd-parity sub-level with higher energy 
and one even-parity sub-level with lower energy. Conse- 
quently, the lowest unoccupied level and the highest oc- 
cupied level have the opposite parities (left panel of Fig. 
2c). As Ad2 is reduced, a level crossing and hence a par- 
ity exchange happens at Ad2 = 2A due to the opposite 
parities of the two crossing levels, leading to a change 
of v from (for the non-coupling evcn-BL films) to 1, 
and hence converting the Z2 number of 2-BL film from 
trivial to nontrivial. The situation, however, is differ- 
ent for the 3-BL film. Now, each 3-fold degenerate level 
is split into three sub-levels, always with one odd-parity 
sub-level sandwiched by two even-parity sub-levels (right 
panel of Fig. 2c) . Consequently, both the lowest unoc- 
cupied and highest occupied levels have the same even 
parity. As Ad2 is reduced, a level anti-crossing instead 
of crossing happens and the number of odd- and even- 
parity sub-levels for both the occupied and unoccupied 
bands remain unchanged, resulting in no change to the 
non-trivial Z2 topology of the 3-BL film. 

In the above simple picture, we have implicitly used 
the prerequisite condition that the level crossing or anti- 
crossing happens only at the F point but not at the M 
point and involves only a few levels close to the band 
gap, which can only be satisfied by an inter-BL cou- 
pling that is not too strong. Also, the coupling can 
not be too weak in order to move away from the weak 
coupling regime that is adiabatically connected to the 
zero coupling limit. Therefore, the intermediate coupling 
strength is a mandatory condition for the non-trivial Z2 
topology. Such " intermediate coupling principle" may be 
utilized in search for the 2D TI phase in other materials. 

As a summary, the film below 4 BLs is an intrinsic 
2D TI with the band structure consisting of " molecular 
orbital" levels without distinction of surface bands from 
bulk bands, as shown in Fig. 3a (left panel) and 3b (left 
panel) for the 2-BL and 3-BL film, respectively. Above 
4 BLs, the band structure is made of surface bands su- 
perimposed onto a 2D projected bulk band, as shown in 
Fig. 2a for the 5-BL film. The projected 2D bulk bands 
keep the non-trivial topology of a 2D TI with a sizable 
gap, but the surface bands gradually appear in the mid- 
dle of the projected bulk band gap with the increasing 
film thickness (see also 15|)) leading to a semiconductor 
(1- to 4-BL) to a semimetal (5- to 8-BL) transition. The 
trivial metallic surface states can be removed by surface 
H adsorption (Fig. 2d) , which effectively converts the Bi 
films into true 2D TIs. 

Our finding of all the ultrathin Bi (111) films being 
2D TIs independent of thickness may provide a possi- 
ble explanation of the recent observation of ID topolog- 
ical metal on the Bi (114) surface [20l|. and will stimu- 
late more experimental interest in this intriguing system. 
The Bi (111) films above 6-BL have already been success- 



fully grown via molecular beam epitaxy 2l|, and hope- 
fully even thinner Bi films can be grown in the near fu- 
ture. The physical mechanism we identified for retaining 
the stable nontrivial Z2 topology has broad implications. 
Most importantly, tuning the inter-layer coupling to the 
intermediate regime, so as to remove the odd-even oscil- 
lation of topological triviality, can be applied as a general 
strategy to obtain the TI phase. 
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